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Abstract 

We study the entanglement between the two beams exiting a Mach-Zehnder 
interferometer fed by a couple of squeezed-coherent states with arbitrary 
squeezing parameter. The quantum correlations at the output are function of 
the internal phase-shift of the interferometer, with the output state ranging 
from a totally disentangled state to a state whose degree of entanglement is 
an increasing function of the input squeezing parameter. A couple of squeezed 
vacuum at the input leads to maximum entangled state at the output. The 
fringes visibilities resulting from measuring the coincidence counting rate or 
the squared difference photocurrent are evaluated and compared each other. 
Homodyne-like detection turns out to be preferable in almost all situations, 
with the exception of the very low signals regime. 
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I. INTRODUCTION 



The notion of entanglement is an essential feature of quantum mechanics, and is strictly 
connected with the nonlocal character of the theory. A two-part physical system prepared 
in an entangled state is described by a non factorizable density matrix. This gives raise 
to partial or total correlation between the outcomes of measurements performed on the 
two parts, even though the parts may be so far apart that no effects resulting from one 
measurement can reach the other part within the light cone. 

Sources of entangled states are required for fundamentals tests of quantum mechanics, 
as well as for applications such as quantum computation and communication [[I]], and tele- 
portation In recent years, entangled photon pairs had been used to test non- locality 

of quantum mechanics 0-0 by Bell inequality [§]. In practice, all the available sources of 
two-mode entangled states are based on the process of spontaneous down-conversion, taking 
place in nonlinear crystals |J. Recently, it has been demonstrated that a beam splitter 
can split an incident photon into two correlated secondary photons [|1^,|11|]. However, such 
process occurs at very low rate, and thus it is of no interest in practical applications. 

In order to study quantum correlations between two radiation modes, and to compare 



different sources of correlated states, one needs to quantify the degree of entanglement |L2 
A good theoretical measure of correlations has been introduced by means of Von-Neumann 
entropy. The entropy of a two-mode state g is defined as 

S[g] = -Tr {glogg} , (1) 

whereas the entropies of the two modes a and b are given by 

S[g a ] = -Tr a {g a \ogg a } S[g b ] = -Ti b { g b log g b } . (2) 

In Eq. (§) g a = Tr b {g} and g b = Tr a {g} denote the state of a and b respectively, as obtained 
by tracing out the other mode from the total density matrix. Following Refs. [|12| -[T5|] we 



define the degree of entanglement of the state g as the normalized excess entropy ]16|] 

e = qr A s[g a ] + S[g b ] - S[g]\ , (3) 

where S[£ TH ] = log(l + N) + AHog(l + iV -1 ), with N = (a^a), denotes the entropy of 
a thermal state, namely the maximally disordered state at fixed intensity. The use of e 
formalizes the idea that the stronger are the correlations in the two-mode state, the more 
disordered should be the two modes taken separately. If g is a pure state, we have that 
S = and S a = S b [IS], so that e = S[g a ]/ S[g aTH ] ranges from zero to unit. Notice that for 



pure state e represents the unique measure of entanglement JL7|] . 

From the experimental point of view, the entanglement can be detected by non-classical 
interference effects occurring in intensity-dependent measurements. In experiments involving 
photon pairs from parametric down-conversion these effects occur when coincident photons 
are mixed at a beam splitter |19|-21|. The probability amplitudes for pairs from the two 



arms show destructive interference, leading to a suppression of the coincidence counting rate 
between detectors surveying the two arms p2| , |23[| . Recently, the spatial effects in two-beam 
interference have been also studied for partially entangled photon pairs f24 |. In the case of 
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more excited states, many photons are present and the connection between entanglement and 



coincidence rate is less transparent. The issue has been received attention P5| , p6| , though a 
general theory has not been developed yet. 

In this paper we study the generation and the detection of entangled states at the output 
of a Mach-Zehnder interferometer fed by a couple of uncorrelated squeezed- coherent states. 
The scheme may be of interest as the output state can be arbitrarily high excited, instead of 
having two photons only. In addition, the degree of entanglement can be tuned by varying 
the degree of squeezing of the input beams, or the internal phase shift of the interferometer. 
As regards the detection scheme, we show that the coincidence counting rate between the 
two output arms corresponds to low fringes visibility. Therefore, we consider instead another 
intensity dependent quantity, namely the squared difference photocurrent, which shows high 
visibility of fringes for the whole range of input squeezing parameter. 

In Section || we study the dynamics of the interferometer, and evaluate analytically 
the degree of entanglement at the output as a function of the squeezing fraction of the 
input beams and the internal phase-shift of the interferometer. In Section |IT| we analyze 
the interference effects occurring in the measurement of the photon coincidence rate and 
of the squared difference photocurrent. The evaluation of the fringes visibility for both 
measurements shows that homodyne-like detection is preferable in almost all situations, 
with the exception of the very low signals regime. Section [1^ closes the paper with some 
concluding remarks. 



II. ENTANGLEMENT AT THE OUTPUT OF A MACH-ZEHNDER 

INTERFEROMETER 

The Mach-Zehnder interferometer we are dealing with is depicted in Fig. |l]a. The input 
signal modes are denoted by a and b, whereas BS\ and BS2 are symmetric beam splitters. 
We also assume that equal and opposite phase-shifts are imposed in each arm of the 
interferometer. The evolution operator of the whole setup can be written as 

where 

U = exp{i^((fb + tfa)} , (5) 

denotes the evolution operator of a symmetric beam splitter. After straightforward algebra 
one rewrites Eq. (|J) as 




which shows that a Mach-Zehnder interferometer is equivalent to a single beam splitter BS^ 
of transmissivity r = cos 2 |, preceded and followed by rotations of ir/2 performed on one of 
the two modes (see Fig. |l]b). 

We consider the Mach-Zehnder interferometer fed by a couple of squeezed-coherent states 

\i> IN )D a (a)D b (a)S a (()S b (()\0) = . (7) 
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In Eq. (|7|) D a (a) = exp(aat — aa) is the displacement operator and S(Q = exp [1/2 (( 2 a^ 2 — 
C 2 a 2 )] is the squeezing operator, |0) denotes the electromagnetic vacuum. There is no need 
to consider a phase-shift between the input modes, as it can be reabsorbed into the internal 
phase shift 0. Without loss of generality, in the following we will consider a complex field 
amplitude a G C and a real squeezing parameter ( e r 6 1. 
The state exiting the interferometer is given by 

IVw) = v MZ {4>) |V%> • (8) 
By exploiting the vacuum invariance Vmz(0)|O) — |0) and using the relation 

exp j— i— cJa j D{a)S{r) exp |z— a^a| = D(—ia)S(—r) , (9) 

we can write |t/> ut) as 

I^out) = exp \i\h%} %b a {a)D h (-ia)Ul U^S a (r)S b (-r)t/l |0> , (10) 

where denotes the evolution operator of the equivalent beam splitter BS^. The "dis- 
placing" part of Eq. fllpp , together with the rotation on the mode b, can be easily rewritten 
as 

exp {i|& f &} U+b a (a)b h (-ia)Ul = b a (ae^)b b (ae^-^) exp {i|& f &} , (11) 

whereas the "squeezing" part needs a little more algebra: specializing a result from Ref. [27 
we can write 

U*S a (r)S b (-r)Ul = expjcos0 ^r(a t2 - a 2 - 6 t2 + b 2 ) + sin0 [r{aW - ab)] J . (12) 

It is worth noting that squeezing at the input is essential to obtain entanglement at the 
output. In fact, for the input state being a couple of coherent states IV'in) = D a (a>)D b (/3)\0) 
the output state is given by I^out) = D a (a cos — ij3 sin (p)D b {— ia sin + a cos0)|O), which 
again is a couple of factorized (uncorrelated) coherent states for any value of the internal 
phase-shift of the interferometer. Actually, the absence of output correlations is due the the 
Poissonian statistics of the coherent states, which implies the absence of intensity fluctuations 



Let us first consider the situation = |. In this case, the transformation in Eq. (|12l) 
reduces to the two-mode squeezing operator S^ 2 '(r) = exp |r(a^ — ab)}, so that the output 
state coincides with a displaced and rotated twin-beam state 



|t/> out > = b a (ae l i)D b (ae^-^)ex P [i^bj |^ TWB ) , (13) 
where the explicit expression of the twin-beam state I^twb) is given by 

1 oo 

IVWb) = S^(r)\0) = -^-^tanhV \k,k) . (14) 

fc=0 
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In order to evaluate the degree of entanglement of |V>out) we use the parameter e introduced 
in Eq. (|3|). The partial trace over a mode, say b, is given by 

1 oo 

Q a = Tr 6 {|Vw)(Vw|} = — rx- Vtanh 2fc r D{ote^)\k){k\b\ae^) , (15) 

cosh r f — ■ 

k=0 

which is diagonal in the basis of displaced number states |^ n ) = D(ae l ^)\n). The set 
of |^ n )'s constitutes an orthogonal basis for the Hilbert space of harmonic oscillator, and 
therefore the entropy S[g a ] can be evaluated as 



S[g a ] = - y^Pn^ogPn p n = tanh 2ri r . (16) 

cosh r 

n=0 

After straightforward calculation we arrive at 



S[g a ] = log(l + v 2 ) + v 2 log(l + -) , (17) 



1 

~ 2 

where v 2 = sinh 2 r is the squeezing energy of each input beam. Notice that S[g a ] in Eq. ([Tj 
is equivalent to the entropy of a thermal state with v 2 photons. The degree of entanglement 
is given by 

c _ log(l +7 jV) +7 jVlog(i + -^) 

log(l + N) + N log(l + i) ' ( } 

where = (a^a) = \a\ 2 + z/ 2 is the total energy of each input signal, and 7 is the squeezing 
fraction, namely the percentage of the total energy engaged in squeezing photons, v 2 = ^yN. 
From Eq. (|T8D it is apparent that the degree of entanglement is an increasing function of the 
squeezing fraction, and that the maximum entangled state (e = 1) at the output is reached 
for a couple of squeezed vacuum (7 = 1) at the input. 

For = 0, the transmissivity of the whole device is equal to unit, and we have I^qut) = \ip m ). 
Therefore, the entanglement is equal to zero, as the input state consists of a couple of 
uncorrelated signals. 

For 7^ 0, 1 it is convenient to evaluate the output state and the entanglement by evolving 
the two-mode Wigner function, which is defined as follows 



W(x a , y a ] %b, Vb) = d/j, a dv a \ d/j, b / du b exp {2i(u a x a - fi a y a + v h x h - fi b yb)} X 
Jr Jr Jr Jr 

x Tr {gDai^a + iv a )D b {n b + iv b )j . (19) 
The ±7r/2 rotations of mode b correspond to simple rotations in the sole b- variables 

g' = e l ^ge~ l % ==> W'(x a , y a ; x b , y b ) = W(x a , y a ; y b , -x b ) 

g' = e"^ g e*' ==► W'(x a , y a ; x b , y b ) = W(x a , y a ; -y b , x b ) , (20) 

whereas the action of the beam splitter BS^, i.e. g' = U^glf^ corresponds to a mixing of 
variables of the two modes, in formula 
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W'(x a , y a ; x b , y b ) = W(x a cos 5 - x b sin 5, y a cos 5 - y b sin 5; 

x a sin 5 + x b cos S, y a sin 5 + y b cos 5) , (21) 



where we use the notation 5 = 4>/2. Using Eqs. ( p20|) and (|21|) the Wigner function at the 
output results 

'f OUtO^oj Dai %bi Db ) = W m (x a cos 5 — y b sin 5, y a cos 5 + x b sin 5; 

cos 5 — y b sin 5, x a sin 5 + y b cos 5) , (22) 

where PUiN^a, y a ; ^6, yb) is a product of two identical single-mode Gaussian Wigner functions, 
corresponding to the couple of input squeezed- coherent states: 

W m (x a , y a ; x b , y b ) = exp i - 2e~ 2r (x a - Re[a]) 2 - 2e 2r (y a - Im[a]) 2 

- 2e- 2r (x b - Re[a]) 2 - 2e 2r (y b - Im[a]) 2 j . (23) 

By the integration over the 6-variables 

''OUT (^tti Da) = 

dx b dy b W 0VT ( 

Veil %bi Vb) , (24) 

Jr Jr. 

and inserting Eqs. (22) and (|23|) in Eq. (24) we obtain 



W (x v) 1 exof ( x -~ Re ^}) 2 {to ~ Im^]) 2 ] 

which represents the Wigner function of the sole mode a after partial trace over the mode 
b. The quantities T, x and T, y in Eq. (|25"D are given by 

S 2 = e 2r cos 2 5 + e- 2 '-sin 2 5 

£ 2 = e- 2r cos 2 <5 + e 2r sin 2 5 , (26) 

whereas is given by 



a^l + ^sin 2 0. (27) 

In order to evaluate entanglement, we note that any unitary transformation T acting on the 
single mode a does not change the value of the entropy |25[| , i.e. S[g a ] = S[Tg a T^}. Using this 



property, we displace with amplitude a^, and then squeeze with parameter r* = log a/S^/S^ 
the Wigner function in Eq. (|2"5"D, thus arriving at the following entropy-equivalent state 

WUx a ,y a ) = exp (-4^1 • (28) 



Remarkably, the Wigner function in Eq. (|28|) coincides with the Wigner function of a thermal 
states with thermal photons given by 
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N+ = \ P y z x - 1] = i 



1 + sin 2 sinh 2 2r 



(29) 



The corresponding entropy can be easily computed, and thus the entanglement at the output 
is given by 

_ log(l + i^) + iVog(l + i) 
6 ~ log(l + iV) + iVlog(l + ^) ' (30) 

As it is expected, one has = for = 0, and = 'yN for = tt/2. 

In Fig. |^a we show the degree of entanglement as a function of the squeezing fraction 7 
and the internal phase-shift 0, in the case of input beams with average photons N = 3 each: 
at fixed 7 the output state ranges from a totally disentangled state for = 0, to a state 
whose degree of entanglement is given by Eq. (|IH) for = |. The degree of entanglement is 
an increasing function of the squeezing fraction 7, with the condition = | corresponding 
to maximum value. Different values of the intensity N does not substantially modify the 
behavior of e versus 7 and 0. In Fig. 0b we report e as a function of the intensity iV 
for different values of the squeezing fraction 7, and for fixed value = | of the internal 
phase-shift: For 7 = 1 one has e = 1 independently on N, whereas for 7 < 1 the degree 
of entanglement becomes a slightly increasing function of N. For highly excited states the 
entanglement is given by the asymptotic formula 

+ (3D 
log N 

So far we have considered the two input states having the same degree of squeezing. However, 
a pair of input states with different squeezing fractions does not substantially modify the 
picture. In this case, in fact, the entanglement still oscillates from e = to a maximum 
value as a function of the internal phase-shift of the interferometer. On the other hand, this 
maximum value is now a function of both the squeezing fractions, and maximally entangled 
states at the output cannot be achieved if one of the input signals is only partially squeezed. 
In Fig. ^ we report the maximum entanglement at the output (obtained for = vr/2) as a 
function of the squeezing fractions 71 of one of the beams for different values of the squeezing 
fraction 72 of the other beam. The plots refers to a situation in which both input beams 
have an average number of photons equal to N = 3. As it is apparent from the plots, the 
output entanglement is an increasing function of both the two squeezing fractions. The 
extreme case in which one of the input signals is no squeezed at all corresponds to a value 
of e always lower than 50%. 



III. ENTANGLEMENT AND FRINGES VISIBILITY 

In this Section, we study the visibility of the interference fringes that are observed, by 
varying the internal phase-shift 0, in intensity measurements at the output of the inter- 
ferometer. In analogy with experiments involving correlated photon pairs, we consider the 
detection of the coincidence counting rate at the output, namely of the fourth-order corre- 
lation function (if) 0l]T \a^a b'b\ip OUT ). However, as we will show in the following, this corre- 
sponds to low fringes visibility, and thus we sought for a more sensitive kind of measurement. 
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The homodyne-like detection of the output difference photocurrent {ip 0VT \a^a — b^b\^ 0VT ) is 



widely used in interferometry P9|-pf|], and generally results in a very sensitive measurement 
scheme. Starting from this consideration, we suggest the squared difference photocurrent 
OUT |(a+a — 6+6) 2 |VWr) as a suitable fourth-order quantity to be measured at the output of 
the interferometer. 

Besides being originated by interference effects, the variations in the quantities mea- 
sured at the output also reflect the variations in the quantum correlations between the two 
output signals. Therefore, the visibility of the interference fringes provides a measure of en- 
tanglement, and comparing the visibility of different measurement schemes provides a way 
to compare their ability in monitoring the variations of quantum correlations between the 
output signals. 

As already mentioned, here we consider the measurement of the coincidence counting 
rate 

K{<j>) = (ip 0VT \a^a 
and of the squared difference photocurrent 
H{4>) = (VootK^o - 

After some algebra, we arrive at the explicit expressions in terms of the input fields 

K(<f>) = VMz(<t>)<Ja tfbV MZ (<f>) = 

= sin 2 5 cos 2 6 [a t2 a 2 + 6+ 2 6 2 + a) 2 b 2 + 6 t2 a 2 ] + (sin 2 6 - cos 2 5) Va 6+6 + 
+ % sin 5 cos 3 5 [ab ]2 b + a t2 a6 - a ] b ] b 2 - a ] a 2 b ] ] + 

+ i sin 3 5 cos 5 [a+a 2 b ] + aW - a ]2 ab - ab t2 b] , (34) 



H{<t>) = Vl iz m^a - tfb) 2 V MZ (4>) = -2K(4>) + 

+ [(a+a) 2 + (6 f 6) 2 ] (sin 4 5 + cos 4 5) - 2 sin 2 6 cos 2 5 [a ]2 b 2 + b ]2 a 2 - a ] a - b ] b] + 
+ 2i sin 5 cos 3 5 [a)a 2 b ] - a ]2 ab + a ] b ] b 2 - a¥ 2 b] + 

+ 2isin 3 5cos<5[a6 t2 &-aW + a t2 a&-a t a 2 & t ] , (35) 



where again we used the notation 8 = ~. Using Eqs. ( p4|) and (^) we are able to evaluate 
the fringes visibility of both detection schemes 

K — K ■ H — H 

t 7- J1 mm J1 min r ±J -max ±± min /o/^\ 

Vk = k Tk — = Tr — IF" • (36) 

-L^-max ~r ^min -^max T -"min 

In Fig. fl we report Vk and Vh as a function of the intensity iV for different values of the 
input squeezing fraction 7. The H-measurement visibility Vh is larger than Vk in almost 
all situations, with the exception of the very low signals regime, where very few photons 
are present. The behavior of fringes visibility versus intensity N also confirms that Vh 
represents a good measure of the entanglement at the output. As it happens for the degree 
of entanglement, in fact, a couple of squeezed vacuum at the input corresponds to maximum 
visibility Vjj — 1 independently on the intensity. On the other hand, the coincidence 
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counting rate shows a visibility Vk that rapidly decreases versus N, and saturates to a value 
well below 1/2. For non unit squeezing fraction, and moderate input intensities iV < 10, the 
behavior of Vjj looks qualitatively similar to that of the degree of entanglement (compare 
Fig. |]b and Fig. |||b), whereas again Vk rapidly decreases. Remarkably, for highly excited 
states N > 10, the visibility Vh has the same asymptotic dependence of the degree of 
entanglement e, in formula 

e^l + fi, (37) 
log N ' V ; 

where the proportionality constant A (7) ~ 1/5 log 7 is roughly proportional to that appear- 
ing in Eq. (H). 



IV. CONCLUSIONS 



The generation and the detection of optical entangled states are important issues, both 
required for fundamentals tests of quantum mechanics, as well as for possible applications. 
In this paper we have studied the entanglement between the two beams exiting a Mach- 
Zehnder interferometer fed by a couple of squeezed-coherent states with arbitrary squeezing 
parameter. The degree of entanglement at the output has been analytically evaluated, as a 
function of the input intensity and squeezing fraction, and of the internal phase-shift of the 
interferometer. Our results indicate that entangled states of arbitrary large intensity can be 
produced by varying the input energy, whereas the degree of entanglement can be tuned by 
varying the input squeezing fraction, and the internal phase-shift. 

An experimental characterization of the output entanglement can be obtained through 
the measurement of the squared difference photocurrent between the output modes. The 
interference fringes that are observed by varying the internal phase-shift show, in fact, 
high visibility for the whole range of input squeezing parameter. 
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FIGURES 




FIG. 1. In (a): schematic diagram of a Mach-Zehnder interferometer. BS± and BS2 are 
symmetric beam splitters, whereas a and b denote the input signal modes. Equal and oppo- 
site phase-shifts are imposed in each arm. In (b): an equivalent scheme for the Mach-Zehnder 
interferometer depicted in (a): a single beam splitter BS^ of transmissivity r = cos 2 J preceded 
and followed by rotations of ^ performed on one of the two modes (here b). 
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(a):the degree of entanglement as a function of the squeezing fraction 7 and the internal 
j>, in the case of input beams with N = 3 each: for fixed 7 the output state ranges from 



FIG. 2. 

phase-shift 

a totally disentangled state for <p = 0, to a state with a degree of entanglement given by Eq. (18) 
for (j> = ^. The degree of entanglement is an increasing function of 7, with </> = § corresponding to 
maximum value. Different values of N do not substantially modifies the picture, (b): the degree of 
entanglement as a function of the total input energy N for different values of the squeezing fraction 
7, and for fixed value <fi = j of the internal phase-shift. 
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FIG. 3. Output entanglement for input signals with different degree of squeezing. The maxi- 
mum entanglement at the output (for <p = ir/2) is reported as a function of the squeezing fractions 
7i of one of the beams for different values of the squeezing fraction 72 of the other beam. Both 
input beams have an average number of photons equal to N = 3. The output entanglement is 
an increasing function of both the two squeezing fractions. The extreme case in which one of the 
input signals is no squeezed at all corresponds to a value of e always lower than 50%. 
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FIG. 4. Fringes visibility as a function of the intensity N for different values of the input squeez- 
ing fraction 7. In (a) the visibility of K-measurement Vk, and in (b) the visibility of H- measurement 
Vh ■ In both plots we report the visibility versus N for five values of the input squeezing fraction. 
From bottom to top we have the curves for 7 = 0.2, 0.4, 0.6, 0.8, and 1.0. As it is apparent, Vh is 
larger than Vk in almost all situations, with the exception of the very low signals regime. 
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